The Polyakov loop and the hadron resonance gas model 



E. Megfacl 

Grup de Fisica Tebrica and IFAE, Departament de Fisica, 
Universitat Autdnoma de Barcelona, Bellaterra E-08193 Barcelona, Spain 

E. Ruiz ArriolcQ and L. L. Salcedc@ 

Departamento de Fisica Atomica, Molecular y Nuclear 
and Instituto Carlos I de Fisica Tedrica y Computacional, 
Universidad de Granada, E-18071 Granada, Spain. 
(Dated: March 2, 2013) 

The Polyakov loop has been used repeatedly as an order parameter in the deconfinement phase transition in 
QCD. We argue that, in the confined phase, its expectation value can be represented in terms of hadronic states, 
similarly to the hadron resonance gas model for the pressure. Specifically, L(T) ~ J Y,a8a e~ Aa / T , where g a 
are the degeneracies and A a are the masses of hadrons with exactly one heavy quark (the mass of the heavy 
quark itself being subtracted). We show that this approximate sum rule gives a fair description of available 
lattice data with Nj = 2 + 1 for temperatures in the range 150MeV < T < 190 MeV with conventional meson 
and baryon states from two different models. For temperatures below 150MeV different lattice results disagree. 
One set of data can be described if exotic hadrons are present in the QCD spectrum while other sets do not 
require such states. 
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Introduction. — The transition from the hadronic phase to 
the quark-gluon plasma phase has been a recurrent topic in 
hadronic physics In pure gluodynamics or, equivalently, 
for infinitely heavy quarks, this is a true phase transition. The 
order parameter is identified as the thermal Wilson line or 
Polyakov loop [2-4], 
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where Aq is the gluon field, T is the temperature and P de- 
notes path ordering. L(T) changes abruptly from zero to 
near N c (the number of colors) due to the breaking of the 
center symmetry Z(N C ) for T > T c , In QCD, i.e., for dy- 
namical quarks, the center symmetry is explicitly broken by 
the quarks and one has instead a smooth crossover Q5|] and 
the critical temperature T c is usually defined by the condition 
L"(T C ) = 0. Lattice simulations show that chiral symmetry is 
restored when quarks and gluons are deconfined. These theo- 
retical insights strongly suggested the experimental quest for 
the quark-gluon plasma in current facilities. The Polyakov 
loop also serves as a gluonic effective degree of freedom 
in the successful Polyakov-Nambu-Jona-Lasinio (PNJL) and 
Polyakov-quark-meson models (PQM) to describe hot and/or 
dense QCD 

Since hadrons (and possibly glueballs) are the physical 
states in the confined phase it should be expected, by quark- 
hadron duality, that physical quantities admit a representation 
in terms of hadronic states. The QCD pressure presents a 
prime example of this, through the hadronic resonance gas 
model (HRGM) {HHH, 
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with g a the degeneracy factor, £ a = ±1 for bosons and 
fermions respectively, and M a the hadron mass. These res- 
onances are the low-lying states listed in the review by the 
Particle Data Group (PDG) ifTsh . Actually, in the large N c 
limit this expectation becomes a true theorem in QCD since 
the flavor resonances become narrow F/M = 0(1 /N c ) (see 
e-g- ll^Eolo . After some controversies, lattice calculations 
seem to suggest that this is also a good approximation in the 
physical case N c = 3 12111 . This problem has been addressed 
within a strong coupling expansion for heavy quarks in 022ll . 

Despite its prominent theoretical role, Q.(x) does not appear 
to be directly accessible in the laboratory, being most naturally 
defined in the im aginary -time formalism of field theory at fi- 
nite temperature 11231 12411 (see Il25ll for its definition within the 
real-time formalism). The realization of the Polyakov loop 
through a static (or heavy) quark and its relation with the 
heavy-quark self-energy (or even with the binding energy be- 
tween a static and a dynamical quark fl26J]), is not new 
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but the phenomenological consequences of this fact have not 
yet been extracted at a quantitative level. Here we argue that 
the Polyakov loop in the confined phase can also be repre- 
sented in terms of hadronic properties in a direct and quanti- 
tative way, similarly to the HRGM for the pressure. 

Polyakov loop and hadronic spectrum. — In the Hamilto- 
nian formulation of QCD lEilEsill . the gauge is partially fixed 
by the condition Aq = and the dynamical degrees of freedom 
are contained in the spatial gluons A and the quarks. The time- 
independent gauge transformations g(x) are still a residual 
symmetry of the Hamiltonian acting in the Hilbert space J$? of 
functionals ^(A^q^q). The gauge group SU(N C ) decomposes 
34? into invariant subspaces labeled by an irrep r at each point 
x: Jif — (B{ r (x)} /-(*)}■ I n tne Euclidean lattice formulation, 
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the role of the integration over Aq, or equivalently integration 
over £l(jt) with the Haar measure, is to project onto the phys- 
ical subspace, which requires a color singlet at every point x. 

An infinitely heavy quark (of a new flavor) sitting at xq is 
a spectator with spin and color degrees of freedom only. For 
the gluons and dynamical (as opposed to heavy) quarks, this is 
equivalent to living in the subspace JK(x Q )=3 : a color singlet 
at every point except xq, which is in the fundamental represen- 
tation (3 for three colors). The projector onto this subspace is 
obtained by adding the factor N c tr(£l(xo)) to the Haar mea- 
sure I130n . For the expectation value of the Polyakov loop this 
immediately implies the relation j3ll l32ll 



ibare(r) = 



nW«-" /r ) 

2Tr phys ( e -"/r) 
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The factor N c in the projector represented the trivial degen- 
eracy of the system formed by gluons plus dynamical quarks 
in the fundamental representation at Jto, and is canceled when 
this is combined with the spectator quark to form a color sin- 
glet. The factor 1/2 removes the double counting from the 
two spin states of the spectator quark. The l.h.s. is indepen- 
dent of the heavy quark spin (as the Polyakov loop carries no 
spin) and this is fully consistent with the well-known heavy- 
quark spin symmetry present in QCD ll33l 34]. The (infinite) 
mass of the spectator quark is not included in H. 

Eq. (0 is exact for the bare Polyakov loop and the partition 
functions in J^hys an d <^rfrc )=3 on me lattice. In the renor- 
malized continuum limit the relation still holds, after remov- 
ing the additional specific UV divergence introduced by the 
heavy quark self-energy in L(T) and Tii UXo (e~ H / T ). Such re- 
moval leaves a nonperturbative ambiguity by an additive con- 
stant in the Polyakov loop free energy F(T) = —TlogL(T) 

The (renormalized) partition functions in Eq. (|3} are sat- 
urated by states of the spectrum. Since the spectator quark 
can be reached smoothly by taking the infinite mass limit of a 
heavy quark at rest, this implies 
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where ra/, denotes the heavy quark mass. Here, the sum in 
the denominator is just the QCD partition function and so 
it includes all possible states made of gluons and dynami- 
cal quarks (labeled by a). On the other hand, the sum in 
the numerator includes all possible QCD states with exactly 
one heavy quark h at restQplus gluons and dynamical quarks 
(jointly labeled by ha). The difference E/ ja — m/, explicitly 
removes the heavy quark mass from the total energy of the 
state. 

For temperatures well below the crossover, we expect the 
previous states to be of hadronic type (and possibly glueballs). 



1 Such statement would be meaningless for dynamical quarks, but not for 
infinitely heavy quarks. 



In particular, the heavy quark h will form a hadron with the dy- 
namical quarks, typically, a meson of hybrid type, i.e., formed 
by the heavy quark and a dynamical antiquark, hq, or a hybrid 
baryon with the heavy quark and two dynamical quarks, hqq. 

The HRGM for the QCD partition function follows natu- 
rally from assuming that the QCD interaction primarily con- 
fines quarks into hadrons and that purely hadronic interactions 
can be neglected. Under this assumption, the numerator of 
Eq. dU contains one hybrid heavy-light hadron at rest plus ex- 
actly the same multi-hadron states occurring in the denomina- 
tor. This yields a cancellation between numerator and denom- 
inator. Therefore, within the same approximations leading to 
the HRGM, we expect the following relation to hold between 
the Polyakov loop expectation value in the confined phase and 
the hadronic spectrum 



L ( T ) ~ 9L Shae 



At a = M ha - m h 



(5) 



Here the sum is over all states made just of a single hybrid 
hadron at rest (with exactly one heavy flavor quark, mass Mj ia 
and degeneracy gi ia ), and no additional hadrons. 

Of course, neither the sum rule in Eq. (0 nor the HRGM 
can be accurate when unconfined states of the spectrum start 
to play a role, that is, for temperatures in the crossover region 
or above. For instance, L(T) is a decreasing function at high 
enough temperatures, in the perturbative regime [39], while 
the Boltzmann distribution in the r.h.s of Eq. is always in- 
creasing as a function of T. On the other hand, the mass of the 
heavy-light hadron is an observable (a renormalization invari- 
ant) but depends on the heavy flavor, while A/ m is universal in 
the heavy quark limit but has some running from ra/,, which is 
itself not an observable. The Polyakov loop is renormalization 
invariant and well defined, modulo the abovementioned shift 
ambiguity in the free energy. This can be compensated by a 
corresponding shift in the heavy quark mass. 

Estimates from the physical spectrum. — In order to com- 
pare different Polyakov loop determinations with the hadronic 
sum rule, they have to be brought to a common renormal- 
ization condition. Two such determinations are related by 
L'(T) = e c l T L{T), for some constant energy shift C. In Fig.Q] 
we compile five Polyakov loop data sets, obtained with physi- 
cal quark masses and three flavors on the lattice IU7M21I1 . The 
plot shows that four of them agree after applying suitable fi- 
nite renormalizations (no attempt has been made to optimize 
the agreement), in a wide range of temperatures. Unfortu- 
nately, the agreement deteriorates at lower temperatures, be- 
low T ps 150MeV. This region is relevant for comparison 
with the hadronic sum rule. Since there is no "true" value 
of L(T), a finite renormalization, or choice of heavy quark 
mass, should be admitted too in the hadronic sum rule. Nev- 
ertheless, large renormalizations (compared with the hadronic 
scale) would be unnatural, and would probably signal an inad- 
equacy of the sum rule, or of the renormalization prescription 
used for the Polyakov loop. A neat way to remove any am- 
biguities is to work with the derivative of riog(L(r)) with 
respect to T. This slope is sensitive to the effective number 
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FIG. 1. (color online). 7Tog(L(r)) as a function of T (units in 
MeV) from simulations on the lattice with 2+1 physical dynam- 
ical quark masses. Lower set of lines: data after a common shift 
C = —150 MeV (just for displaying purposes) for continuum extrap- 
olated stout (black solid line, "circle") 12 ill . HISQ/tree action N t — 12 
scale set r\ (blue solid line, "square") and fa (blue dashed line, "up 
triangle"), and asqtad action N t = 12 scale set r\ (red solid line, 
"down triangle") and fa (red dashed line, "rhombus") 11711. Upper 
set of lines: same data (except asqtad scale set r\ ) with shifts C = 0, 
0, -lOMeV and -30MeV, for stout, HISQ/tree scale r\, HISQ/tree 
scale fa, and asqtad scale fa, respectively. 
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FIG. 2. (color online). Extraction of the slope ^(Tlog(L(r))) (di- 
mensionless) as a function of T (in MeV) from the same four lattice 
data sets of Fig. [T] The error bars in the slope were obtained from 
assuming a linear interpolation between measured points. The black 
thick line indicates the average in the range of temperatures common 
to the four sets. The yellow strip indicates the uncertainty. 



of states at a given temperature. Although with some noise, 
Fig.[2]shows that a signal can be extracted in this way. 

A natural step is to check to what extend the hadronic sum 
rule is fulfilled by experimental states compiled in the PDG. 
Several sources of error should be kept in mind when doing 
this, among other, that not all needed states may have been 



FIG. 3. (color online). Comparison of ^(Tlog(L(r))) (yel- 
low strip) with -jj{T\og{\YlaSha e ^^'""^ T )) from hadronic states 
(mesons plus baryons): lowest-l ying hadrons from PDG (solid brown 
line, label a), RQM states from 14(1 kill with quark c (solid red line, 
label b), and with quark b (dashed red line, label c), and bag model 
estimate including states up to A = 5500 MeV (solid blue line, label 
d). 



compiled, that the heavy quarks in nature have a finite mass, 
that their current mass is scale dependent, and that the quark 
masses on the lattice may not be identical to the physical ones. 
Hadrons with a bottom quark would be optimal, due to the 
large quark mass compared to Aqcd, but the available data 
are scarce, so we turn to charmed hadrons. Specifically, we 
consider the lowest lying single-charmed mesons and baryons 
with u, d, and s as the dynamical flavors, with quarks in rel- 
ative s-wave inside the hadron. For mesons, these are usu- 
ally identified with the states (spin-isospin multiplets) D, D s , 
D* (2010) andDJ, and for baryons, with A e , £ e (2455), E c , S' c , 
Q. c , E c (2520), S c (2645), and H(2770). A total of 12 meson 
states and 42 baryon states. 

The plot in Fig. [3] shows that the lowest-lying states fall 
short to saturate the sum rule, regardless of the choice of mass 
of the charmed quark, m c . This is not surprising as any model 
predicts many excited states on top of the lowest-lying ones, 
as is also the case for light-quark hadrons. Adding more states 
from the PDG does not seem practical due to the fragmen- 
tary information available. Instead we turn to hadronic mod- 
els. The aim is not so much to have a detailed description of 
the various states but to give a sufficiently good overall de- 
scription of the whole spectrum. To this end we consider the 
relativized quark model (RQM) (4(J 41], and the bag model 
142, 43ll . We have verified that the RQM provides a good 
account of the trace anomaly in 112 ill . The total number of 
hadron states computed in 1 4(J 41 ] with one c quark is 1 17 for 
mesons and 1470 for baryons, corresponding to a maximum 
value of A = M — m c about 1500 MeV. For hadrons with one 
b quark, 87 mesonic and 1740 baryonic states, with a similar 
upper bound for A. In both cases we have supplemented miss- 
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ing states with strange quarks by means of the equal spacing 
rule B44I1 and a s quark mass of 109MeV (extracted from the 
lowest-lying hadrons masses). The prediction based on these 
hadronic states is displayed in Fig. [3] The two sum rules are 
closer to the Polyakov loop result but still tend to stay below 
it in the range T < 175 Me V, a consequence of the truncation 
of states to A < 1500MeV. It is noteworthy that the bottom 
sum rule gives a better value, as it would be expected due to 
the larger mass of the b quark. 

The other model we consider is more schematic but allows 
us to easily include a larger number of states. This is a simpli- 
fied MIT bag model 11421 14311 . in order to correctly count the 
number of states without fine details such as multiplet split- 
tings. As bag energy we take 



A = 



Ej/IjCDf — Z 47T 

R ^ 



# 3 B + £m/, 



(6) 



where n„ o, and m; are the occupation number, bag frequency 
and current quark mass (m u = m c / = 0, m s = 109 Me V). The 
model gives directly A after adjustment of R, without center of 
mass corrections, since the heavy quark has actually infinite 
mass (not included), sits at the center and plays no active role. 
The sum over i runs over just one antiquark for hq mesons and 
two quarks for hqq baryons. We set B = (166MeV) 4 j43ll . 
and find m c = 1390 Me V and Z = 0.5 from a fit to the single- 
charmed lowest-lying mesons and baryons masses. The re- 
sult for the bag model with A < 5500 Me V is displayed in 
Fig. [3] The more complete description of the hadronic spec- 
trum gives a better account of the Polyakov loop data in the 
range 145 MeV < T < 175 MeV. We have checked that: i) For 
this range of temperatures heavier states become irrelevant for 
the sum rule. We see this by projecting the cumulative num- 
ber of states assuming a power law dependence for mesons 
and for baryon [ 4511 . ii) Truncation to A < 1500 MeV gives a 
result quite consistent with that of RQM. And iii) A similar 
power-law projection of the spectrum for the RQM, to esti- 
mate the effect of adding the states above A = 1500 MeV, also 
reproduces quantitatively the bag model result. Fig. [3] also 
shows that the hadronic sum rule eventually overshoots the 
Polyakov loop result, as the crossover to the deconfined phase 
sets in. This mimics the same behavior in the HRGM j2lll46ll . 

In Fig. [4] we display lattice data for T\ogL(T) vs the 
hadronic sum rule. The hadronic estimate (with no additional 
finite renormalization in the case of the bag model) describes 
well the various lattice data sets in the range 150 MeV < 
T < 190 MeV, and even lower temperature data from IU7II . 
However, the steeper slope displayed by the stout action data 
112 111 for T < 150 MeV cannot possibly be saturated with con- 
ventional mesons and hadrons, since all these states have al- 
ready been accounted for. Inclusion of exotic hadrons, hqqq 
(tetraquarks) and hqqqq (pentaquarks), do actually produce 
the required slope (see Fig. 0J, and it would imply a much 
shorter temperature range of applicability of the hadronic sum 
rule. At present, the various sets of lattice data disagree at the 
lowest temperatures, and also it is unsettled whether exotic 
hadrons are present in the QCD spectrum or not [47]. Our 
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FIG. 4. (color online). T log(L(T)) (in MeV) for the four lattice data 
sets of Fig.Q]compared to the hadronic sum rule from the bag model. 
Solid blue line: estimate from conventional mesons and baryons (C = 
0). Dashed green line: estimate when exotic hadrons are included, 
applying a shift C = -40 MeV. 



analysis implies that resolution of one of these issues would 
shed light on the other. 

Final comments. — From QCD considerations, we derive a 
Boltzmann distribution formula in terms of hadrons for the 
expectation value of the Polyakov loop in the confined phase, 
as required by quark-hadron duality, where its real and posi- 
tive character becomes manifest. Our derivation exposes the 
obvious fact that the Polyakov loop gets its expectation value 
from the dressing with dynamical quarks or antiquarks. Since 
U(Nf) is an exact global symmetry, the numerator in Eq. (|3), 
and hence L(T) itself, can be decomposed into separate con- 
tributions from different flavors and different baryon numbers. 
Such decompositions are in principle accessible to lattice cal- 
culations (although with a difficulty similar to that of introduc- 
ing a chemical potential) and they would provide further in- 
formation about the interplay between the QCD thermal state 
and the heavy quark spectrum. 
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